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Abstract 

We present results for a complete set of polarization observables for jet production in lepton proton 
collision, where the final state lepton is not observed. The calculations are carried out in collinear 
factorization at the level of Born diagrams. For all the observables we also provide numerical 
estimates for typical kinematics of a potential future Electron Ion Collider. On the basis of this 
numerical study, the prospects for the transverse single target spin asymmetry are particularly 
promising. This observable is given by a certain quark-gluon correlation function, which has a 
direct relation to the transverse momentum dependent Sivers parton distribution. 

1 Introduction 

In lepton nucleon scattering one normally detects the scattered lepton in order to determine the vir- 
tuality Q 2 of the exchanged gauge boson. Provided that Q 2 is sufficiently large, one can use, for a 
number of final states, the machinery of QCD factorization (see Ref. [TJ for an overview) in order to 
separate the short-distance physics from the non-perturbative long-distance physics encoded in differ- 
ent parton correlation functions. In data analyses for lepton nucleon scattering, QCD factorization 
not only has been used to get a handle on objects like ordinary forward parton distributions, but also 
to address generalized parton distributions and transverse momentum dependent parton distributions. 

In the present work we study inclusive jet production in lepton proton scattering with the scattered 
lepton going unobserved, i.e., Ip — > jetX. In this case, the transverse momentum of the jet can serve 
as the large scale, which is needed for justifying a calculation in perturbative QCD. The kinematics 
of this process is rather simple — in particular, simpler than the kinematics of semi- inclusive deep- 
inelastic scattering (DIS) — and in essence coincides with the one of, e.g., single jet production in 
hadronic collisions. We focus here on jet production, as opposed to hadron production, because this 
process, in principle, can provide more direct information about parton correlation functions of the 
proton as no uncertainties from parton fragmentation are involved. 

We consider all possible polarization observables for the process Ip —> jetX at the level of Born 
diagrams. Neglecting parity violating effects as well as transverse polarization of the initial state 
lepton, one can identify three spin-dependent cross sections: oil, &ut, and olt-, where the first index 
refers to the lepton polarization and the second one to the proton polarization. While o~ll, like the 
unpolarized cross section auu, is a twist- 2 observable, the latter two are twist-3 effects. 



In order to compute gjjt and ctlt we make use of collinear twist-3 factorization, which was pi- 
oneered already in the early 80's [21 [3] and applied for the first time to single spin asymmetries 
(SSAs) in [2]. These early treatments of twist-3 SSAs were later revisited and improved [HE]. I n the 
meantime many papers on SSAs and related observables in the collinear twist-3 formalism exist; see, 
e.g., Refs. [6l426|. In particular, various works are dealing with the QCD-evolution of the relevant 
twist-3 correlators (see |27H31j and references therein). Moreover, the relation between the collinear 
twist-3 factorization and a factorization in terms of transverse momentum dependent parton correla- 
tors [32Tf34] has been studied in detail for semi-inclusive deep-inelastic scattering and the Drell-Yan 
process [35H38] . 

Based on our numerical estimate for a typical kinematics of the currently discussed/planned Elec- 
tron Ion Collider (see, e.g., Ref. pSJIIO]), the transverse SSA Ajjt = (Jut/ouu appears to be rather 
promising. The QCD-description of Ajjt contains a specific twist-3 quark-gluon-quark correlator - 
the so-called ETQS (Efremov-Teryaev-Qiu-Sterman) matrix element [2J[4]. As pointed out in [4"ljl4"2]. 
the ETQS matrix element is related to the transverse momentum dependent Sivers function [J3]. (For 
experimental studies of the Sivers effect in semi-inclusive DIS we refer to |44}I45]. while extractions 
of the Sivers function from data were discussed in [4*6H5Tj ) . Thus, measuring Ajjt in lp — > jetX 
would give a direct (complementary) handle on the Sivers effect. Our prediction for the longitudi- 
nal double spin asymmetry All = oll/&uu is at the percent level, while for Alt-, computed in a 
Wandzura-Wilczek-type approximation, we obtain only a tiny effect. 

Note also that a detailed study of Aut ialp—^ jet X, based on factorization in terms of transverse 
momentum dependent correlators, can be found in Ref. |52| . That work represents an extension and 
update of a related earlier investigation of Ajjt for Ip — > ir X [53J. Moreover, Ajjt for pion production 
was also considered in the collinear twist-3 approach in two conference proceedings |54j . 



2 Factorization and predictive power 

In this section, we present a QCD factorization formalism for the inclusive high transverse momentum 
jet production in lepton hadron collision and provide brief arguments why this factorization formalism 
should be valid. We also provide the prescription for calculating the short-distance hard parts of the 
formalism, and discuss the predictive power of perturbative calculations. 



2.1 Factorization formulas 

With a large momentum transfer, the transverse momentum Pjt = \Pjt\ of the inclusive jet in 
high energy collisions, the short-distance dynamics takes place at a time scale of 1/Pjt, which is 
much shorter than the typical time scale of hadronic physics, 0(fm). The quantum interference 
taking place between these two very different time scales is likely suppressed by the ratio of these 
two scales. Perturbative QCD factorization of hadronic cross sections effectively neglects the power 
suppressed quantum interference and factorizes the cross section into a product or a convolution of 
two probabilities: one for finding active parton(s) inside the identified hadron(s), and the other is the 
short-distance part of partonic cross section(s). For example, the leading power contribution to single 
inclusive jet production at large transverse momentum in hadron hadron collisions, h{P) + h'(P') — )• 
jet(Pj) + X, can be factorized as [l], 

da hh'^jet(Pj)X /■ /■ d ~ab^jct(Pj)X 

dP JT dy J J dX '^ ^''^ — dP JT dy — ^' X ' ' PjT ' V ' ^ ' ^ 

where 2~^ a b runs over all parton flavors, f^ h (x,fj.) is the (unpolarized) parton distribution function 
(PDF) of flavor a and momentum fraction x of hadron h, and \x is the factorization scale. Since 
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Pjt is the only large observed momentum, Eq. ([TJ is a collinear factorization formalism [TJ. The 

da dPjJdy^ ( x ' x '' PjTi V, A*) i n Eq. (P) is the perturbatively calculable short-distance hard part, which 
is effectively equal to the partonic jet cross section for the collision between two partons a and b with 
all collinear sensitive contribution removed. The predictive power of Eq. ([TJ relies on our ability to 
calculate the partonic hard parts systematically in perturbative QCD order-by-order in a s , and the 
universality of PDFs. The factorization formalism in Eq. ([TJ works extremely well for describing single 
inclusive jet data at the Tevatron for over ten orders of magnitude in the production rate |55j . 

Our ability to calculate the short-distance partonic hard parts in Eq. ([TJ relies on the fact that the 
factorization of short-distance dynamics is not sensitive to the long-distance details of the colliding 
hadron(s). That is, the factorization formalism in Eq. ([TJ, which is valid for two colliding hadrons, 
should also be valid for the collision of two asymptotic partons. By applying Eq. ([1]) to the collision 
of two partons of various flavors, we can derive all short-distance hard parts order-by-order in powers 
of a s |56| . Since the factorization formalism in Eq. ([TJ is not sensitive to the details of the colliding 
particles, we expect that the same factorization formalism in Eq. ([TJ is also valid for single inclusive 
jet production at high Pjt in lepton-hadron collision, 1(1) + h(p) — > jet(Pj) + X, as 

da lh^jct(Pj)X , , d frab-+Jet(Pj)X 

dPjTdy ^gy^A &V)]**h (x, M ) dPjTdy (x,z',PjT,y,ri, (2) 

where J2 a runs over the lepton, the photon, and all parton flavors, while J2b runs over all parton 
flavors, fi^(x,fj,) is the nonperturbative distribution to find a lepton, photon or parton inside the 
colliding lepton with the momentum fraction x. Like in the case of hadronic collisions, the short- 
distance hard part da dPjTdy j) (x, x', Pjt, U, y) can be perturbatively calculated order-by-order in the 
coupling constant by applying the factorized formalism in Eq. ([2]) to the collision between various 
lepton, photon or parton states. For the leading contribution, it might be reasonable to keep the cross 
section at the lowest power in a em while including radiative corrections from the strong interaction in 
powers of a s . 

Like all perturbative QCD factorization approaches, the predictive power of Eq. ([2]) relies on the 
infrared safety of the hard parts and the universality of the long-distance distributions. Unlike the 
hadronic case in Eq. ([TJ, the jet production in lepton hadron collisions requires a set of new non- 
perturbative distributions, /" (x, y) with a = l,~/,q,q,g. The operator definition of the lepton PDFs 
should be the same as the proton PDFs' except that the proton state is replaced by the state of the 
lepton [57]. The operator definition of the lepton distribution inside a lepton is very similar to that 
of a quark distribution, 

f[ /l {x^) = J ^ e ^ + r (Z)Sl |^' (0 )l^ W7 (o, r )/(e-)|z, s i), (3) 

where W 7 (0,£~) = exp[— ie Jq dy~ A+(y~)] is the gauge link in an Abelian gauge for the lepton 
moving in the u +z" direction. If we neglect the role of the strong interaction, we can calculate the 
lepton distribution perturbatively in QED. At the lowest order, f \ ^ l (x,fi) = S l ' l 6(l — x). The operator 
definition of the photon distribution inside a lepton is the same as the operator definition for the gluon 
distribution inside a proton with the color dependence removed and the proton state replaced by the 
lepton state, 

where is the electromagnetic field strength tensor, which can be expressed in terms of gauge 
invariant electric and magnetic fields. The photon distribution of the lepton, fj (x, /i), could be 
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Figure 1: Diagram (a): parton model representation for Ip — > jetX. The jet is produced by the struck 
quark, and the final state lepton goes unobserved. Diagram (b): contribution from quark-gluon-quark 
correlation. This diagram, together with its Hermitian conjugate which is not displayed, needs to be 
taken into account when computing twist-3 observables. 



calculated perturbatively in QED if we neglect the strong interaction. But, in general, it is a non- 
perturbative distribution. 

Based on the same argument that the factorization of the short-distance dynamics is not sensitive 
to the details of the colliding particles, and the fact that the single inclusive jet cross section at 
high transverse momentum has only one observed large momentum transfer Pjt, we expect that the 
twist-3 collinear factorization formalism for calculating the transverse-spin dependent cross section 
in hadronic collisions [2,4,5j can be applied to the single inclusive jet cross section in lepton hadron 
collisions. Therefore, 



dP JT dy 



J dxff(x^) J dx'T h J h {x',x'^,S T )H ah ^ p ^ x {x,x',P JT ,y^), (5) 



where Act(St) = [a (St) — c(— 5r)]/2 is the transverse-spin dependent cross section, with the spin 
vector St = \S~t\ of the transversely polarized colliding hadron. The Tp h (x,x, fi, St) with b = q,q,g 
in Eq. ([5]) is the universal twist-3 parton correlation function relevant for the SSAs [HET], an d the 
jjab^jet(Pj)x ( X)X ' ^ p JT ^y^ ^ with a = l,j,q,q,g and b = q,q,g are the process-dependent short- 
distance hard parts whose leading order contributions are derived below. 

2.2 Partonic hard parts 

In this subsection, we provide the prescription for calculating the partonic hard parts of the factor- 
ization formulas in Eqs. ([2]) and ([5]). 

To calculate the short-distance hard parts, ^—jpj^j^— (x, x', Pjt, y, /•*) in Eq. ((2j), we apply the 
factorization formalism to the collision between all possible combinations of two asymptotic incoming 
lepton, photon, or parton states. For the leading order (LO) contribution, the jet cross section is given 
by the lowest order lepton-quark scattering, as sketched in Fig. QJa), and the jet is effectively given 
by the final-state quark: jet(Pj) —¥ q(Pj) at the lowest order. The corresponding hard part can be 
uniquely derived by applying Eq. ([2]) to the collision of the lepton on a quark state: I — > I and h — > q, 

lq^q(Pj)^jct(Pj) _ Ml ,q/q l q ^ q (Pj)^j e t(Pj) _ ~lq^q(Pj)^jct(Pj) 

"(2,0) ~ J 1(0) 59 J 1(0) 59 "(2,0) ~ "(2,0) ' W 

where the superscript "(2,0)" indicates two powers of a em and zeroth order in a s , and (8) represents the 
convolution over lepton or parton momentum fraction as shown in Eq. ([2]). At the lowest order, the 




(a) (b) 

Figure 2: (a) Leading order tree diagrams for photon scattering on an asymptotic quark state, (b) 
Leading order tree diagrams for photon scattering on an asymptotic quark-gluon composite state. Note 
that only the gluon interaction with the observed final-state parton is nonzero, while the interaction 
with the unobserved parton cancels. 



short-distance hard part for the jet cross section is effectively the same as the lepton-quark scattering 
cross section for producing the quark q(Pj) and is perturbatively finite. 

If we apply the factorization formalism in Eq. ([2]) to photon-quark collision by letting I — > 7 and 
h — > q, as sketched in Fig. [2]^a), we can have two additional tree-level contributions to the jet cross 
section, 



W^q(Pj)^}et(Pj) 



a 



(1,1) 



4?' 
^ 41 



4V 



(7 



a 



jq^q(Pj)->-}et(Pj) 
(1,1) 

~/q^9(Pj)-+}et(Pj) 
(1,1) 



jq-> q (Pj)-+j e t(Pj) 
"(1,1) 
jq-+9(Pj)-+}et(Pj) 
(1,1) 



a 



(7) 



The difference of these two contributions is whether the jet is generated by an energetic quark or 
a gluon. Similarly, if we apply Eq. to photon-gluon collision, we obtain two more tree-level 
contributions to the jet cross section: J ^~~* anc } g.J9^>q{P.j)->}ct(Pj) ^ gj nce fa e photon 



distribution of the lepton f^ l (x,fi) carries at least one power of a em higher than the leading term of 
fi\x,^L) : these photon-parton contributions could be formally considered as an order a s correction 



to the LO term in Eq. 



However, since the photon distribution of the lepton f^ l (x,fi) has a 
large QED logarithm perturbatively, for certain kinematics these terms could be more important than 
typical higher order corrections. 

If we apply the factorization formalism in Eq. ([2]) to parton-parton collision by letting / — > a 
and h — > b with a,b = q, q, g, the partonic hard parts are to be the same as those in the hadronic 
collisions [HE]. Since the parton distribution functions of the lepton fi (x, /x) have at least two powers 
of o, em perturbatively, even the leading Born contribution here should be considered as higher order 
correction in powers of a s . 

Higher order corrections to the short-distance hard parts of the factorization formalism can be 
derived in the same way by applying Eq. ([2]) to the scattering of various partonic states at higher 
orders in a s . For example, we can calculate the next-to- leading order (NLO) contribution a"? 1 



by applying Eq. ([2]) to the lepton-quark collision at order a s 



(2,1) 



a 



Iq-fjet(Pj) 
(2,1) 



Ml 
J 1(0) 



h(o) 



a 



,fj/l gv f q/q 
+/l(l) » Ji(o) 



Iq^jet(Pj) 
(2,1) 



+ fl 



l/l 
(0) 



fill 
/l(l) 



^q^q(Pj) ~iq^g{Pj) 
"(1,1) ^ (1,1) 



lq-+jet(Pj) 
'(2,0) 



(8) 



which can be written as 



~lq^]et(Pj) 
(2,1) 



a 



a 



Iq^jet(Pj) 
(2,1) 



fill 
J 1(1) 



a 



Iq^iet(Pj) 
(2,0) 



fill 
/l(l) 



^q-^q(Pj) , --yq^giPj) 
"(1,1) ^ (1,1) 



(9) 
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Figure 3: (a) Next-to-leading order tree diagrams for lepton scattering on an asymptotic quark state, 
(b) Next-to-leading order tree diagrams for lepton scattering on an asymptotic quark-gluon composite 
state. 



where the first term on the right-hand-side (RHS) is the partonic cross section given by the real 
Feynman diagrams sketched in Fig. EJa) plus the virtual diagrams from one-loop corrections to the 
tree-diagram in Fig. QJa). The real diagrams in Fig. [3ja) have three potential collinear logarithmic 
divergences. One comes from the final-state gluon radiation when the gluon is parallel to the parent 
quark. Such final-state collinear divergence is taken care of by the jet definition and its finite cone size. 
The other two come from the initial-state: when the gluon line is almost parallel to the incoming quark 
or when the photon is about parallel to the incoming lepton. As a result of QCD factorization, these 
two divergences are systematically removed by the second and the third terms in Eq. ([9]) , respectively. 
The hard 2 — > 2 scattering cross section of the second term, £(2~o) ^ is derived in Eq. (El), while 



a 



jq^q(Pj) 
(1,1) 



and a 



(i,i) 



of the third term are given in Eq. (|7|). 



Similarly, we can derive another short-distance contribution at the same order, a l %~$ , by 



applying Eq. ([2]) to the collision between a lepton and a gluon, 



(7 



(2,1) 



Zg-s-jet(Pj) 
(2,1) 



fi/g 
h{i) 



a 



lq->jet(Pj) 
(2,0) 



4) ® 



-■yg^q(Pj) ±rtg->q(Pj) 
"(1,1) ^ (1,1) 



(10) 



where the second and the third terms on the RHS again remove the collinear divergence of the partonic 
scattering cross section. In general, the perturbatively calculated hard parts are effectively equal to 
the partonic cross sections with all collinear divergences removed. 

In the collinear factorization approach, the spin-dependent cross section to the SSAs comes from 
the interference of the real part of the scattering amplitude with one active parton and the imaginary 
part of the scattering amplitude with two active partons [21S1E1E] - For the LO contribution, the 
partonic hard part H lq ^ q ^ Pj ^ x (x, x' , Pjt, y) m Eq. ([5]) is given by the diagram in Fig. DJb), and will 
be calculated in the following section. 

For the higher order corrections, the partonic hard part H ab ^ > ' Jet ^ Pj ^ x (x, x' , Pjt, U, I 1 ) in Eq. ([5]) 
can be calculated in the same way by applying the formalism to the various partonic states. For 
example, the LO photon-parton scattering contribution to the SSA comes from the interference of 
Feynman diagrams in Fig. EJa) and[2](b). Similarly, the NLO lepton-quark scattering contribution to 
the SSA comes from the interference of Feynman diagrams in Fig. El(a) and[3^b). In the rest of this 
paper, we will present our derivation and results of the LO contribution to various spin asymmetries. 
We will leave the explicit treatment of higher order corrections to future work. 



3 Kinematics and analytical results 

In this section we present some details of the kinematics for the process 1(1) + p(P) — > jet(-Pj) + X, as 
well as the tree level formulas for the various observables. We use the momenta of the particles to fix a 
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coordinate system according to e z = P/\P\ = — e x = Pjt/\Pjt\, and e y = e z x e x . Mandelstam 
variables are defined by 

s = (l + P) 2 , t = (P-Pjf, u = {l-Pj)\ (11) 

while on the partonic level one has 

s = (I + k) 2 = xs, t = {k-Pjf = xt, u = (l-Pj) 2 = u , (12) 

with k denoting the momentum of the active quark in the proton; see also Fig. []Ja). The momentum 
fraction x specifies the plus-momentum of the quark through k + = Using s + i + u = one 

finds that x = —u/(s + t). In other words, the longitudinal momentum of the struck quark is fixed 
by the external kinematics of the process, like it is in fully inclusive DIS. Of course, this no longer 
applies once higher order corrections are taken into account. For the numerical estimates we will use 
Pjt, and the Feynman variable xf (defined in the lepton-proton cm- frame) for which one has 

2P.Jz t-u 

XF = f=r = • ( 13 ) 

Vs s 

Next, we turn to the polarization observables for Ip — > jet X, which we compute in the collinear 
factorization framework. We restrict ourselves to one-photon exchange between the leptonic and the 
hadronic part of the process. Allowing for longitudinal polarization of the initial state lepton, as well 
as longitudinal and transverse polarization of the proton target, one finds the following expression for 
the cross section!! 

73^ „,2 a 2 



■2nM e%S l T P 3 JT 



+ \ x 2M S T ■ Pjt 



Tf(x,x) -x4-7f(x,x) 
ax 



j7 Huu 
tu 



~g a (x) - x^~9 a (x)) | H LL + xg a T {x) I 



(14) 



In Eq. (|14p . which is the main analytical result of our work, A; and X p represent the helicity of the 
lepton and the proton, respectively. One can project out the four independent components of the 



cross section in (TH|) according 


to 






&UU = 


\ (<K+, +) + *(-, +) + ") + CT (-> ") 


)■ 


(15) 


&LL = 


±([a(+,+)-a(- +)]-[a(+,-)-a(- - 


-)]), 


(16) 


(jut = 


| (K+, tv) + ty)] " H+, iy) + °{- 




) , (17) 




l(W+Ax) - <r(-Ax)] ~ [<r(+,U) ~ <*{- 


i,)} 


) ■ (18) 


In these formulas, '+' and '— 


' indicate particle helicities, whereas 'tx/y 


( \-x 


iy) denotes transverse 



polarization of the proton along e x / y (—e x / y ). 

As already mentioned, both ajju and otl are twist-2 observables. We computed them on the basis 
of diagram (a) in Fig. [1] by applying the collinear approximation to the momentum k of the active 



J For a generic 4- vector v, we define light-cone coordinates according to i; ± = (v° ± v :i )/v / 2 and v T = (v 1 ^ 2 ). 
2 Polarization degrees are suppressed in the cross section formula (|14[) . 
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quark. In the case of auu the result contains the unpolarized quark distribution /f , while for an the 
quark helicity distribution gf shows up. The hard scattering coefficients for these two terms in (Q~4j), 
expressed through the partonic Mandelstam variables in (|12p . read 

The cross sections ctut an d &lt (3rd and 4th term on the r.h.s. of (|14p . respectively) represent 
twist-3 observables. (Note that M is the proton mass, and e T = e ^ with e 0123 = 1.) The transverse 
SSA Aut = (Tut /cruij is analogous to the SSA An which has been extensively studied in one-particle 
inclusive production for hadron-hadron collisions; see also Ref. [52], and [58-60] for experimental 
results from RHIC. A similar observable was proposed in [6HIS2] for semi-inclusive DIS, but in this 
case the final state lepton still needs to be observed. 

Calculational details for such twist-3 observables in collinear factorization can be found in various 
papers; see, e.g., Refs. [HElEldS]. We merely mention that one has to expand the hard scattering 
contributions around vanishing transverse parton momenta. While for twist- 2 effects only the leading 
term of that expansion matters, in the case of twist-3 the second term is also relevant. In addition, 
the contribution from quark-gluon-quark correlations, as displayed in diagram (b) in Fig. CO needs to 
be taken into consideration. The sum of all the terms can be written in a color gauge invariant form, 
which provides a consistency check of the calculation. 

The quark-gluon-quark correlator showing up in ajjT is the aforementioned ETQS matrix element 
Tp(x, x) [210]. The peculiar feature of this object is the vanishing gluon momentum — that's why it 
is also called "soft gluon pole matrix element" . If the gluon momentum becomes soft one can hit the 
pole of a quark propagator in the partonic scattering process, providing an imaginary part (nontrivial 
phase) which, quite generally, can lead to single spin effects 011]. Note also that in our lowest order 
calculation of ajjT no so-called soft fermion pole contribution (see [20] and references therein) emerges. 
For &LT another quark-gluon-quark matrix element — denoted as g a ; see, in particular, Refs. [8,22,28] 
— appears, together with the familiar twist-3 quark-quark correlator g T . 

We use the common definitions for g\ and gT- The quark-gluon-quark correlators Tp and g are 
specified according tc|l 

T F (x,x) = ^|^^e^-(P,5 r |^(0) 7 + ^ +l (r)^(r)|P,ST)(ie?^) , (20) 

~ g{x) = J_ / ^l e -p+r ( St 

yy ' 2M J 2vr V T 

x <P, <? r |VK0) 7 57 + (iD T -ig J™ dCF+\C]\^C)\P,S T )i (21) 

with F^ v representing the gluon field strength tensor, and = — igA 11 the covariant derivative. 
Equations ()20|) and (|2ip hold in the light-cone gauge A + = 0, while in a general gauge Wilson lines 
need to be inserted between the field operators. 

It is important that Tp and g are related to moments of transverse momentum dependent parton 
distributions. To be explicit, one has [22j l 



7TTp(x,x) = -Jd 2 k T ^fMx,k 2 T )\ Dis , (22) 
9(?) = J d 2 k T ^gMx,k 2 T ), (23) 



3 Note that in the literature different conventions for Tf exist. 
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where we use the conventions of Refs. [63H65] for the transverse momentum dependent correlators f^p 
and g\T- In Eq. (14) we take into account that the Sivers function f^j, [43] depends on the process 
in which it is probed |66|J67]. In order to obtain numerical estimates for ajjT and ult we will exploit 
the relations in (122]), (US). 



4 Numerical estimates 

Now we move on to discuss numerical estimates for the polarization observables. To this end we 
consider the three spin asymmetries All, Ajjt, and Alt, whose definitions are repeated here for 
convenience, 

A LL = ^, A UT = ^, A LT = ^. (24) 
<?uu &UU &UU 

To compute cjjju we use the unpolarized parton distributions from the CTEQ5-parameterization [68] . 
The helicity distributions entering oll are taken from the GRSV-parameterization [69]. For the 
ETQS matrix element Tp we explore two choices: (1) we use the relation (|22[) between Tp and 
the Sivers function, and take from the recent fit provided in Ref. [51]; (2) we use Tp from the 
parameterization obtained in [U] — taking into consideration the recently discovered sign change [70] 
— by fitting transverse SSAs measured in hadronic collisions. 

In the case of <jlt one needs input for gp and g. For gx we resort to the frequently used Wandzura- 
Wilczek approximation [71] (see [72] for a recent study of the quality of this approximation) 

g T (x)* t-gx{y), (25) 
Jx y 

whereas for g we use (|23p and a Wandzura-Wilczek-type approximation for the particular A^-moment 
of gix in ((23]) [73] , leading to 

n dy 



g(x) ~ x / —gi(y). (26) 
Jx y 

We mention that (|26p and a corresponding relation between chiral-odd parton distributions were used 
in |74|,I75| in order to estimate certain spin asymmetries in semi-inclusive DIS. The comparison to data 
discussed in [75] looks promising, though more experimental information is needed for a thorough test 
of approximate relations like the one in (|26p . Measuring the double spin asymmetry Alt, which we 
consider in the present paper, may provide such a test. 

For our numerical estimates we use leading order parton distributions, and take into account 
the three light quark flavors. The transverse momentum of the jet Pjt serves as the scale for the 
parton distributions. For the following reasons the scale-dependence of all the asymmetries is rather 
weak: because the leading order evolution kernels for fi and g\ are the same, All is almost scale- 
independent. This also applies to Alt when using the approximations ([25]) and ([26]) . Since both the 
parameterization of the Sivers function in |51j as well as the one for the ETQS matrix element Tp 
in [9"][70] are related to the unpolarized distribution f\, also the scale dependence of Ajjt is quite mild. 

Our results are for typical kinematics accessible at a potential future Electron Ion Collider |39|.I40| : 
we consider the energies ^fs = 50GeV and yfs = 100 GeV. The asymmetries are either presented as 
function of xp for fixed Pjt or vice versa. 

We start by discussing the twist-2 asymmetry All- As shown in Fig. [4] this observable is relatively 
small (on the percent level). It is largest in the backward region (negative xp), and rises with increasing 
Pjt- Despite the small effect, measuring All could provide complementary information on the quark 
helicity distributions of the proton. On the other hand, for the longitudinal double spin asymmetry 
in Ip — > jetX one faces the same problems one has in inclusive DIS: quarks and antiquarks enter with 
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_j 0.03 r 




Figure 4: All as a function of xp (left) and Pjt (right). Solid line: yfs = 50GeV; dashed line: 
ys = 100 GeV. 
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Figure 5: Ajjt as a function of xp for y/s = 50 GeV (left) and y/s = 100 GeV (right). The solid line is 
for Tp taken from the fit of the Sivers function in [51] and using the relation in (|22p . while the dashed 
line is computed with Tp from [9jl7"0]. 



equal weight, and a flavor separation is hardly possible. However, if instead one considers An, for 
inclusive hadron production these problems can, in principle, be circumvented like in semi-inclusive 
DIS. According to Fig. [4j All clearly increases towards lower values of yfs. Therefore, All for 
Ip — > H X (at y/s < 50 GeV) should definitely be a very interesting observable for studying the quark 
helicity structure of the proton. 

Let us now turn to the transverse SSA Ajjt, which is displayed in Fig. [5] and Fig. [6j For both 
parameterizations we obtain a healthy asymmetry in the forward region, with effects at the level 
5 — 10%. (Note also that, for the parameterization taken from [51], our numerical results in the 
collinear approach are similar to those obtained in Ref. [52] by using factorization in terms of transverse 
momentum dependent parton correlators and the same input for the Sivers function.) The asymmetry 
drops with increasing Pjt and hardly changes when varying y/s (see Fig. [6]). The weak energy- 
dependence of Ajjt appears because, according to (fT4|) . both ajju and ajjT have the same hard 
scattering coefficient. From Fig. [5] it is obvious that the transverse SSA Ajjt seems very promising 
in order to experimentally constrain the ETQS matrix element Tp and the Sivers function. Such a 
constraint is of utmost importance as the existing parameterizations even differ in sign [70J. If Pjt 
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0.08 




Figure 6: Aut as a function of Pjt at xf = 0.1. Solid line: \fs = 50 GeV; dashed line: \fs = 100 GeV. 
The matrix element Tp is taken from the fit of the Sivers function in [51] and using the relation in (1221) . 




Figure 7: Alt as a function of xf (left) and Pjt (right). Solid line: \fs = 50 GeV; dashed line: 
y/s = 100 GeV. Note that Alt is scaled by a factor of 100. 



is sufficiently large, our leading order calculation should give a reliable estimate of the asymmetry. 
One may even be able to check whether Tf(x,x) changes sign as function of x, as has been recently 
speculated [76] , 

Finally, our numerical estimates for Alt are shown in Fig. This asymmetry is apparently 
too small to be measured. The hard scattering coefficient of glt is the same than the one for oll- 
Therefore, Alt decreases with increasing energy like All does. The main reason for Alt being even 
much smaller than All can be traced back to the factor x showing up on the RHS of the Wandzura- 
Wilczek-type approximation ([26]) . Because of this, an experimental study of Alt can serve as an 
interesting check of the relation (|26p . as already pointed out above. In addition, like in the case of 
All-, measurable effects for Alt hi hadron production at lower values of yfs can be expected. 

5 Summary and discussion 

We have studied a complete set of polarization observables for the process Ip — > jetX. Neglecting 
parity violating contributions and transverse polarization of the lepton one can consider three spin 
asymmetries: All, Aut-, and Alt- We have computed these asymmetries at the level of Born diagrams 
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in collinear factorization. Moreover, numerical estimates for typical kinematics of a potential future 
Electron Ion Collider have been provided. (We have explored the cm-energies i/s = 50GeV and 
\fs = 100 GeV.) In the following we summarize our findings and add some discussion: 

• We have discussed in detail how to calculate the process systematically at higher orders in 
perturbation theory. Mainly for two reasons it is important to extend our explicit calculations 
to the 1-loop level. First, one must investigate how stable the various asymmetries are upon 
inclusion of NLO-corrections, in particular, the logarithmically enhanced 7 — q channel as pointed 
out in Sec. [2j Second, a 1-loop calculation for twist-3 observables, in which derivatives of quark- 
gluon-quark correlators (see dTp/dx and dg/dx in (fbi|) ) appear, has never been done before. 

• The numerical result for the double spin asymmetry All (Alt) is small (tiny) — on the percent 
level for All. However, in both cases significant effects can be expected for lower values of y/s, 
let's say around 10 — 20 GeV. Of course, in this region one cannot perform jet-measurements 
but has rather to consider inclusive hadron production. 

• We find that in the forward region the transverse SSA Ajjt can become of the order 5— 10 %. This 
observable gives a direct handle on the ETQS twist-3 matrix element Tp — from a theoretical 
point of view it is one of the simplest observables for addressing Tp (to leading order neither 
soft fermion poles, nor tri-gluon correlations, nor hard gluon poles contribute) — and also, by 
means of (|22p . to the transverse momentum dependent Sivers function fyj,. Given the fact that 
at present even the sign of Tp is unclear [70] , experimental information on this observable would 
be very valuable. 

In general, we believe that there is sufficient justification for further exploring the potential of high- 
energy lepton nucleon scattering with an unidentified final state lepton. Such type of reaction may 
also constitute an interesting part of the physics program at a future Electron Ion Collider. 
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